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Abstract. This paper deals with Mellin transforms of rational functions g/f 
in several variables. We prove that the polar set of such a Mellin transform 
consists of finitely many families of parallel hypcrplancs, with all planes in 
each such family being integral translates of a specific facial hyperplane of 
the Newton polytope of the denominator /. The Mellin transform is naturally 
related to the so called coamoeba A'j := Arg (Zf), where Zf is the zero locus of 
/ and Arg denotes the mapping that takes each coordinate to its argument. In 
fact, each connected component of the complement of the coamoeba A'^ gives 
rise to a different Mellin transform. The dependence of the Mellin transform 
on the coefficients of /, and the relation to the theory of A-hypergeometric 
functions is also discussed in the paper. 



1. Introduction 

The Mellin transform Mh of a locally integrable function h on the positive real axis 
is defined by the formula 

/*°° r] 

(1) M fc (a)= / h(z)z s -, 

Jo z 

provided the integral converges. Here s is a complex variable s — a+it. The Mellin 
transform is closely related to the Fourier-Laplace transform via an exponential 
change of variables. More precisely, the value of M^(s) is equal to the Fourier- 
Laplace transform of the function x i-> h(e~ x ) evaluated at the point —is. 

In this paper we consider Mellin transforms of rational functions h = g/f, 
where g and / are polynomials. Since the general case is easily settled once we 
have fully investigated the special case where g = 1, and since this will simplify 
our notation and therefore clarify our argument, we shall focus mainly on the case 
h = l/f. 

Let us start by considering the one-variable situation. Given a polynomial 

f(z) = a + aiz + . . . + a m z m 

we assume for the moment that its coefficients a ,...,a m are positive numbers. 
Then the integral with h = l/f converges and defines an analytic function in 
the vertical strip < a < m. 

One can in fact make a meromorphic continuation of this Mellin transform and 
write it as 



(2) 



M 1/f (s) = $(s)T( s )r(m-s), 
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where $ is an entire function. To see this, let us first look at the case of a simple 
fraction l/(a + 6a;). In this case one has the explicit formula 

M 1/f ( S )=a s - 1 b- s r( S )T(l- S ), 

which can be easily established for instance by means of a residue computation. 
Now, considering a general product 

m 

one can decompose 1// into a sum of simple fractions and hence immediately 
deduce that its Mellin transform will be of the form \&(s)r(s)r(l — s), for some 
entire function VP. In fact, a straightforward residue calculation shows that 

*( s ) = - e -™^res[^- 1 //(^)], 

and by the theorem on the total sum of residues it then follows that 

*(1) = #(2) = . . . = *(m - 1) = . 

This means that we obtain formula (J2J) with $(s) = ^(s)/[(l — s) • • • (to — 1 — s)]. 

We have thus found that all the poles of the meromorphic continuation are 
located at the two integer sequences 0, — 1, —2, . . . and m+ 1, m + 2, . . . emanating 
from the end points of the interval [0, m]. Notice that this interval is the Newton 
polytope of our one-variable polynomial /. 

As a matter of fact, in the above discussion we did not actually need to assume 
that the coefficients ao, 01, . . . , a m be positive. A necessary and sufficient condition 
for the argument to work, and in particular for the integral to converge, is that 
f(z) 7^ for all real positive values of z. Another way of formulating this latter 
condition is that / should have no roots with argument zero. 

We now turn to the multidimensional case, and we begin looking at a simple 
example with the denominator / being an affine linear polynomial. 

Example 1. Consider the polynomial f(z) = 1 + z\ + z 2 . The Mellin transform of 
the corresponding rational function 1// is then given by the integral 

poo poo si S2 J- J„ poo /.oo pea 

/ / Zl Z2 dZldZ2 = / / / z s S- 1 z?- 1 e-< 1+ *i+^dz 1 dz 2 dv, 

Jo Jo 1 + z l + z 2 z x z 2 J Q J Jo 
which after the coordinate change t — vz\, u = vz 2 becomes 

/•OO />OG /'OO 

/ t^e-'dt u S2 - x e- u du v- Sl - S2 e- v dv = r(s 1 )r{ S2 )T(l-s 1 -s 2 ). 
Jo Jo Jo 

We shall see in this paper that the fact that the poles of the Mellin transform 
are determined by a product of T-functions is not unique for the special cases we 
have considered so far. In fact, the Mellin transform of a rational function in any 
number of variables will turn out to be always a product of T-functions in linear 
arguments, multiplied by some entire function, so that the only poles of the Mellin 
transform are the poles of the T-functions. Moreover, the configuration of polar 
hyperplanes is governed by the Newton polytope of the denominator polynomial, 
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with one family of parallel hyperplanes emanating from each facet of the Newton 
polytope. Many of the results have been announced previously in [13] . Let us note 
in passing that a similar phenomenon can be observed also for Mellin transforms 
of more general meromorphic functions, with transcendental denominators. As an 
illustration of this we recall the classical formulas 

M 1/e * = r(s) and M 1/{e *_ x) = ((s)T(s) = (s - l)((s)T(s - 1) , 

where £ denotes the Riemann zeta function. 

2. Newton polytopes and (co)amoebas 

Throughout this paper / will denote a complex Laurent polynomial 

(3) f(z) = ^2 a a z a , a„eC«, 

where A C Z n is a finite subset and C* denotes the punctured complex plane C\{0}. 
Here we use the standard notation z a — z" 1 • • • z" n for z = (zi, . . . , z n ) £ C™. 

The Newton polytope A f of the polynomial / is defined to be the convex hull 
of A in R n . We shall primarily be interested in the case where Af has a nonempty 
interior. Like any other polytope, the Newton polytope Af may be alternatively 
viewed as the intersection of a finite number of halfspaces: 

v 

(4) A f = f){<je« n ; (Mfc,«r> >^fc}, 

k=l 

where the fit £ 2™ are primitive integer vectors in the inward normal direction of 
the facets of Af, and the Uk £ Z are integers. 

In general we will let T denote a face of the Newton polytope of arbitrary 
dimension, < dim(r) < dim(A/), and we define the relative interior relint(r) of 
such a face to be the interior of T viewed as a subset of the lowest dimensional 
hyperplane containing it. For each face F we also introduce the corresponding 
truncated polynomial 

fr = ^2 aa7 ' a ' 

consisting of those monomials from the original polynomial / whose exponents are 
contained in the face F of the Newton polytope Af. 

The amoeba Af and the coamoeba A'f of a polynomial / are defined to be the 
images of the zero set Zf = {z £ C" ; f(z) = } under the real and imaginary 
parts, Log and Arg respectively, of the coordinatewise complex logarithm mapping. 
More precisely, one has 

A f = Log(Zf) and A) = Arg(Z / ) , 

where Log(z) = (log \zi\, . . . , log \z n \) and Arg(z) = (arg(zi), . . . ,arg(z„)). 

Writing w = x + iO € C™ and z = Exp(w) = (exp(uii), . . . , exp(w n )), one 
obtains the identities x — Log(z) = Re (w) and 9 = Arg(z) = Im (w), as illustrated 
in the following picture. 
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(M/27rZ) n 

Figure 1 . Real and imaginary parts of the complex logarithm mapping 

The amoeba Af is a subset in W 1 , whereas the coamoeba A'* can be viewed as 
being located cither in the n-dimensional torus (M/27rZ)" or as a multiply periodic 
subset of R n . This reflects the multivaluedness of the argument mapping. 

For brevity of notation we denote the amoeba and the coamoeba of a truncated 
polynomial /p by Ar and A' T . 




3. Mellin transforms of rational functions 

The natural generalization to several variables of the standard Mellin transform of 
a rational function 1// is given by the integral 

f z s dzi A ... A dz n f e< a,a! > , 

(5) Mi//(s) = 4 rn zi-z* =Lw) dxiA -- A x - 

where K™ = (0,oo)" denotes the positive orthant in K". In order for such an 
integral to converge one has to make some assumptions about the exponent vector 
s and also about the denominator /. It turns out that it is not enough to demand 
only that / be non- vanishing on K™ . 

Definition 1. A polynomial f is said to be completely non- vanishing on a set X if 
for all faces T of the Newton polytope A / the truncated polynomial fr has no zeros 
on X . In particular, the polynomial f itself does not vanish on X . 

Remark. This concept of completely non-vanishing polynomials is closely related 
to the notion of quasielliptic polynomials discussed in [5]. 

Theorem 1. If the polynomial f is completely non-vanishing on the positive orthant 
K™ then the integral converges and defines an analytic function in the tube 
domain {s£ C™ ; Rc s = cr G int A f } . 

Proof. It will suffice to prove that for any given s with a g int Af there are 
positive constants c, k > such that 

(6) ^(e^e-^l = \f{e x )\e~ {a ' x) >ce k ^, i£l". 

The proof is by induction on the dimension n. The case n = 1 is easy. Let a 
and fi with a < (3 be the two endpoints of A f . Then for sufficiently large negative 
x one has 

|/(e x )|e— > l -\a a \e^~ a ^\. 
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and for sufficiently large positive x 

\f{e x )\e-°- x > i^leC-^*!. 

Now make the induction hypothesis that the inequality Q holds for dimensions 
< n — 1, and consider a polynomial / of n variables. For each face V of Af, with 
< dimT < n — 1, the given point a can be expressed as a convex combination 

a = A or + (1 — A)r r , 

where err S relint(r) and rr € rclint(conv(A \ T)). Fix a choice of such a point ar 
in each face T, and consider for each T the new convex polytope 

A r = conv((A \ T) U cr r ) ■ 

Notice that when dimT = 0, that is, when T is a vertex of Af, one has Ar = A/. 
Notice also that the original point a belongs to each Ap- 

Let Cr be the outer normal cone to Ar with vertex at err: 

(7) C r = {xeR n | (Z-tTr,z-(rr) <0, V£e A r } . 

All these cones Cr are of full dimension n and together they almost cover the entire 
space M. n . More precisely, the complement 

R n \(\jc r 

r 

is a bounded subset of E™ . Then one can let Cr be a slightly smaller closed convex 
cone, still with vertex at ar, such that Cr \ or is contained in the interior of Cr, 
and such that the complement of the union UrCr is still a bounded set. Notice 
that for x £ Cr \ or the inequality in ^ will be strict, and we may in fact assume 
this to be true uniformly. 

We now observe that it is enough to prove the estimate ^ for x € Cr- Actually, 
it suffices to do it for ieCr\ -Bfl(O) for some large ball £?i?(0). From the induction 
hypothesis we conclude that there are constants cr such that 



fr(e*)e-^,x) 



> c r > , x € 



Indeed, fr{e x ) is a function depending on fewer variables than n, since it is homo- 
geneous in directions orthogonal to T, and ar € relint(A/ r ). 

For each face T let gr{z) be the function containing all the monomials not on 
T so that fr + gr = /• Now we use the decomposition / = f r + gr so that one 
obtains 

(8) f{e x )e- {a ' x) = e^"^ {fr{e x )e- {aT ' x) + gr(e x )e~( ar ^) . 

Take x e C r and write x = a T + y. Recall that a € A/. The first factor e ^r-<y,x) 
can be estimated from below by Coe k ' y ' with the positive constants cq and k given 
by Co = exp(oT — a, ar), and 

k = min{(cr r - a,y) ; \y\ = 1 , er r + y G Cr } > . 

Assuming, which we may, that > |ot|, & n d hence that |x| — |ot| > 1^' — o"r| = V, 
we find 

e {a r -<r,x) > ^gfc^l ^ 

where ci = Coe _fe ' <Tr '. 
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To finish the proof of the inequality we now only need to bound the expression 
in brackets in (JsJ) from below by a positive constant. From the induction hypothesis 
we have that 

\fre- {a ^ x) \ >c r > 0, 
and it is therefore enough to show that the remainder term gr{e x )e~(' 7l '' x ' > stays 
small, say < cp/2. We have the identity 

9r(e x )= E a « e{a ' X) E ~ a «z {a ' v) ■ 
aeA\r aeA\r 

Since a € Ar we have a strictly positive constant 

k a = min{ (or - a, y) ; \y\ = 1 , or + y € C r } , 

and hence 

\a a e {a ' x) | = |o a e (CTr - Q ^ | < |o Q | e"*" M . 
This means that for some large enough R one has 

\gr{e x ) e~ {ar ' x) \ < c r /2 , whenever |or + x\ > Rq , x e C r . 

Hence there is an inequality l/Xe 31 ) exp (— (or, sc))| > cr/2, and we can conclude 
that for all x in Cr \ -Bh(O), for some large ball -Br(O), one has the desired estimate 



f(e x ] 



with c = ciCr/2. □ 

Having thus established the convergence of the integral ^ defining the Mellin 
transform, we now turn to the question of finding its analytic continuation as a 
meromorphic function of s in the whole complex space C™. The polar locus of 
the meromorphic continuation turns out to be a finite union of families of parallel 
hyperplanes. The normal directions of these hyperplanes are precisely the vectors 
fj,k from the representation Q of the Newton polytope A f . 

Theorem 2. If the polynomial f is completely non-vanishing on the positive orthant 
R™ and its Newton polytope Aj is of full dimension, then the Mellin transform Mil t 
admits a meromorphic continuation of the form 

N 

(9) M v/ ( a ) = $(a) nr((w,s)-n), 

k=l 

where $ is an entire function, and where Hki v k are the same as in equation Q). 

Before giving the proof of this theorem let us illustrate the idea of the argument 
by means of a specific example. 

Example 2. Consider the polynomial f(z) = 1 + Z2 + z\ + z\z\. It is easy to check 
that the representation ^ of its Newton polytope is given by 

Wi > 0} n {o-i - o- 2 > -1} n {-2^ - ct 2 > -4} n {o- 2 > 0} , 

so in this case the Newton polygon A f has four inward normal vectors given by 
^1 = (1,0), fi 2 = (1,-1), fta = (-2,-1), and ^ = (0,1). 

We know from Theorem 1 that the Mellin transform M 1 /j is holomorphic for all 
s = (si,s 2 ) whose real part cr = (01,02) lies inside the Newton polygon Af. In 
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order to achieve a meromorphic continuation of M%/f across the left vertical edge 
of Af it suffices to perform an integration by parts with respect to Z\ . Indeed, this 
gives us the identity 



(10) 



M 1/f (s) = 



(2zf + zizf) z s x z^ dz\dz2 
Jo (! + z 2 + z\ + z x z 2 ) 2 z x z 2 



and we claim that this integral, that is, the Mellin transform multiplied by s x , 
converges for all s with real part a in the dark triangle on the left in Figure 2. 
This means that M x /f has been continued meromorphically over the hyperplane 
si = as desired. To verify the claim we decompose the integral in (10) into two 
Mellin type integrals containing the integrands 2z 2+Sl z s 2 2 / f 2 and z| +Sl z 2+S2 / f 2 
respectively. Since the Newton polygon of the denominator f 2 is equal to the 
original Af dilated by a factor 2, we see that the convergence domains for these 
two integrals are given by the translated polygons (—2, 0) + 2A/ and (—1, — 2) + 2A/ 
respectively. The sum of the integrals therefore converges on the intersection of the 
translated polygons, and this is precisely the dark triangle on the left in Figure 2. 







/ \ 






\ 


1 \ 






\ 




Figure 2. The convergence domains (dark) of the integrals after 
the two cases of integration by parts, given as the intersection of 
two translated copies of Af2 = 2A/. The dashed polygon is Af. 



We have thus seen how a meromorphic continuation can be carried out in the 
horizontal direction, that is, in the direction given by /ii. Suppose next that we 
wish to obtain a similar mermorphic extension across the upper left edge of Af, 
the one with normal vector /12 = (1,-1). The way to acheive such a "directional 
integration by parts" is to suitably introduce a parameter A and then to differentate 
with respect to A. More precisely, we make the coordinate change z\ 1— > \z\, 
z 2 1— > \~ 1 Z2 and obtain 

11 ioio A + z 2 + \ 3 zj + z x z 2 z x z 2 

Here the left hand side is obviously independent of A. Hence so is the right hand 
side, and after differentiating and plugging in A = 1 we find that 



= (1 + s x - 8 2 ) 



dz x dz 2 /• 00 / >0 ° (1 + 3zf) z{ x z s 2 2 dz x dz 2 



f(z) z x z2 Jo Jo f{z) 2 z x z 2 
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This relation can be re-written as 



M 1/f (s) = r 



(1 + Zzl) z{ x z s 2 2 dzxdz-i 



S2 Jo JO 



/(*) 



z\z 2 



and reasoning as above we find that this latter integral converges for all s with real 
part a in the dark polygon on the right in Figure 2, thereby yielding a meromorphic 
continuation across the hyperplane Si — S2 = — 1. 

This method of repeatedly performing integration by parts in all the directions 
Hk, by using the corresponding coordinate changes Zj X^ kj zj, is the basis for our 
proof of Theorem 2, and it gives a global meromorphic continuation of the original 
Mellin integral. For our special example, the picture below indicates the full set of 
polar hyperplanes, going out in all directions from the Newton polytope A/. 




Figure 3. Polar hyperplanes of the Mellin transform M\/f. 



Remark. For the Mellin transform of a general rational function gj f each mono- 
mial in the numerator g produces an integral similar to the one in the theorem, 
except that we get a shift in the variable s by an integer vector. This corresponds 
to a translation of the Newton polytope of /, and hence also of the domain of 
convergence of that particular integral. If g has several monomials it can very well 
happen that the intersection of all the corresponding shifted polytopes is empty. In 
that case the integral defining the Mellin tranform may not actually converge for 
any values of s. Nevertheless, performing the meromorphic continuation of each of 
the integrals associated with the monomials from g and then summing these mero- 
morphic functions, we still obtain a natural interpretation of the Mellin transform 
Mg/f as a meromorphic fucntion in the entire s-space. 

PROOF. We prove that the integral ^ can be re-written in such a way as to 
make it have a larger convergence domain, at the expense of having to multiply the 
integral with reciprocals of linear terms corresponding to the poles of the gamma 
functions. In order to achieve this we shall repeatedly "integrate by parts" in each 
of the directions given by the vectors /z/.. Each such step consists in first making 
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the corresponding dilation {z\, . . . , z n ) H> (A Pfcl zi, . . . , X^ kn z n ) of the coordinates, 
then differentiating with respect to the dilation parameter A, and finally setting A 
equal to 1. 

Note that, if T is the facet of with inward normal vector fi}-, the truncated 
polynomial /r has the homogeneity fr(X^ k z) = X Vk fr{z). Hence, the scaled poly- 
nomial X~ Uk /(A^*" z) has the property that all its monomials with exponents from 
r have coefficients that are independent of the parameter A. This means that in 
the differentiated polynomial 

there are no monomials with exponents from the facet T. Its Newton polytope is 
therefore strictly smaller than A f , with the integer from the original inequality 
(/Ltfe,<r) > Uk being replaced by Vk + 1, or possibly by an even larger integer. 

Starting from the original integral expression (J5j) for the Mellin transform Mj / f , 
introducing the parameter A, and keeping in mmd that Miff itself is of course 
independent of A, we obtain 

= 4r [ [ X " kZr ^ = 4- f " dZ 



dX f{X^ k z) z dXl J R n X- Vk f{X^"z) z 

which upon performing the differentation and setting A = 1 yields the identity 

(11) ({Vk,s)-v k )M 1/f (s)= f 

As we shall iterate this procedure it will be important to keep track of polytopes 
of different sizes, and to this end we introduce, for any vector 7 € Z™, the notation 

N 

a( 7 ) = f| e Rn ; ^7*}- 

k=l 

In particular, we have Af = A(i/). Now let m £ N N be a given vector, and perform 
the integration by parts rrij times in the direction of fij, for each j = 1, . . . ,N. The 
total number of such integrations will thus be \m\ = mi + ... + 771^. We claim that 
this iterative process leads to an expression for the Mellin transform that is of the 
form 

1 f z s g m (z) dz 



<12) M "' [a) = W^A-)lnf^' 

where g m is a polynomial whose Newton polytope satisfies A 9m C A(|m|i/ + m) 
and Uj(s) = Jl"=o s ) — v j + w ith the convention Uj = 1 if nij = 0. 

The proof of the claim is by induction. First we check that it holds true in the 
case \m\ — 1, that is, when m is a standard unit vector with 1 in the A:'th entry 



and zeros elsewhere. Indeed, this is precisely the content of formula ( 11 ), where we 
recall that the Newton polytope of g ek is contained in A(v + e^). 
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Assume now the claim to be true for some given vector to, and let us show that 
it then holds also for m! = m .+ ej., where e k is a unit vector as before. Introducing 



again the dilated coordinates X^ k z, we can re- write the integral in equation (12 1 as 



z s \-\™Wk-mkg m (\»k z } dz 
A-(i+MK/(A^z) 1+ l m l V 



We should then differentiate this expression with respect to A and put A = 1. 
When the derivative falls on the monomial in front of the integral we get a factor 
(Hk , s) — Vk + rrik which is precisely what needs to be incorporated into the function 
itfc, and when we differentiate under the sign of integration we arrive at an expression 
of the form 

z s g m '(z) dz 
/(z) 2 +M z ■ 



The new polynomial in the numerator is g rn i (z) = (l + \m\)g ek (z)g m (z) — j '(z)g m {z) , 
where 



~g m (z) = —(\-\ m ^- m *g m {\^z)) 



To finish the proof of the claim we must show that A 9m , C A(|to'|^ + to'). We shall 
use the fact that the Newton polytope of a product of two polynomials is equal 
to the (Minkowski) sum of their Newton polytopes, and also the obvious general 
inclusion A(7) + A(S) C A(7 + 5). Recalling the induction hypothesis, we first 
see that the Newton polytope of the product g ek gm is contained in the polytope 
A(i/+efc)+A(|m|z/+m) C A((l + \m\)v+m+e k ) = A(\m'\v+m'). Then, since the 
polynomial g m has no monomials with exponents on the plane (/x^., a) — |m|z/fc+rrifc, 
we similarly get that the Newton polytope of the other term fg m is contained in 
A(^) + A(|m|z/ + m+ e^.) C A(|m'|^ + m'). From this the claim follows, that is, 
the Mellin transform is given by (12 1 with g m satisfying A 9m C A(|m|^ + m). 



Our next step is to prove that the integral in ( 12 ) converges and defines an 
analytic function for all s with real parts a in the enlarged polytope A(y — m). 
By considering separately each term of g m , we can infer from Theorem 1 that the 
domain of convergence will contain (the interior of) the intersection 



(13) 



n Ed 



\m\)A 



f 



r] 



of translates of dilated copies of A/ . 

Let us check that A{y — m) is indeed a subset of (13). Take an arbitrary 
<tq G A{v — m). By definition it satisfies the inequalities 



(14) 



k = 1, 



,iV. 



and 



In order to see that cto also belongs to the intersection (13), take any r £ A 9n 
observe that the polytope (l + |m|)A/ — ris given by the inequalities 

(15) (/*k,ff + T) > (l + |m|)i/ fc , k = l,...,N. 

What we have to show is that cto satisfies these inequalities. In view of the inclusion 
A Bm C A(|m|z/ + to), we have (/ife,T) > \m\vk + for all k. Together with (14) 
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this gives 

(/ifc,a- + t) = (nk,a ) + (/Xfe,r) > v k - m k + \m\v k + m k = (1 + \m\)v k , 



the intersection ( 13 ) 



so <tq does indeed satisfy (151, and since r was arbitrary it follows that erg li es hi 



In the interior of the domain A(v — m) + i R™ the only poles of M x /f are given 
by Uj{s) = 0, j = 1, . . . , iV. All these poles are simple. This is the same polar locus 
as for the product FT. T((/j, k , s) — v k ). By the theorem on removable singularities it 
follows that the quotient M 1 /f/\\ k r({/j, k ,s) — v k ) = $ is holomorphic for er inside 
the polytope A(j/ — m). But here m £ N N is arbitrary, and since the union of all 
the A(V — m) is the entire space K™, we conclude that $ is in fact an entire function 
as claimed in the theorem. □ 



4. Two special cases 

In certain situations we are able to make our description of the Mellin transform 
even more precise, and explicitly compute the entire function $ that occurs in 
front of the gamma factors in Theorem 2. We have already encountered such a 
case in Example 1 of the introduction, where we considered the transform of the 
simple fraction 1/(1 + z\ + 22)- Elaborating this example just a little further, and 
considering a more general linear fraction l/(co + C\Z\ + ... + c n z n ) with each 
coefficient c k being a positive real number, one easily deduces the formula 

(16) M l/f { s ) = <# + ~+'*- 1 c[ si ■ ■ ■ c- s - r( Sl ) • • • r( s „)r(i - 8l - . . . - s n ) . 

So in this case the entire function $ is equal to the elementary exponential function 
s 1 — ^ Cq 1+ " +s " _1 c^ Sl ■ • • c~ Sn and in particular different from zero everywhere. 

We shall now consider two families of examples that both generalize the case 
of a linear fraction, namely products of linear fractions and rational functions that 
are obtained from linear fractions by means of a monomial change of variables. 

Proposition 1. Assume that the polynomial f(z) = YiT^oi^ ( afc ' z )) * s a V r °duct 
of affine linear factors, with each a k £ MS_. Then the Mellin transform of the rational 
function l/f is equal to 

(17) M 1/f ( s ) = $(«)r(«i) . . . r(« n )r(m + 1 - *i -...-*„) , 

with the entire function $ given by 

dri ■ ■ ■ dr m 



Here a m denotes the standard m-simplex { r 6 IR™ ; T k < 1 } , and the ckfc(r) are 
affine linear forms defined by 

(ai(r), . . . ,a„(r)) = (l - ^V fe )a + nai + . . . + T m a m . 

PROOF. We begin by first computing the Mellin transform of a power of the 
type 1/(1 + (c, z)) m+1 . By performing repeated integrations under the sign of 
integration we get 



(-1)™ d m 



A=l 
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Then, recalling the formula ( 16 1 and using the simple identity 



(-1)™ d m _ ASl+ ... +s „_i _ 1 T(m + 1 - s n ) 



A=l ml r(l — Si — ... — s n ) 



we find that 

M 1 /( 1+ ( C)a ))m+ l = cp 1 • • • c~ s "T(si) • • • T(s n )T(m + 1 - si - . . . - s n ) 
Next we make use of the generalized partial fractions decomposition 

1 /" dn ■ ■ ■ dr„ 



which occurs in the theory of analytic functionals and Fantappie transforms, see 
for instance pQ or . From this formula we immediately obtain 



which yields ( 17). □ 



In particular, when m = n = 1 and f(z) = (1 + aoz)(l + a±z) we obtain the 
entire function 

*(«) = /\ i °l~ a - a l s ^J_y Tes[z s^i /f{z)] 

Jo ((1 - r)a + rai) 1 - s a x - a 1 - s ' 

in accordance with the formulas mentioned in the introduction above. Similarly, 
when n = 2 and /(#i, Z2) = (1 + aoiZi + 002-22) (1 + Bn^i + ^12^2) the entire function 
becomes 

$ , . = r 1 rfr 

io ((1 - r)a i + Tan) Sl ((1 - r)a 02 + ra^)* 2 
Here one may remark a close connection to the classical Euler beta function B. 
Namely, if we let the coefficients a^2 and an become zero, we are left with 

$0i,s 2 ) = aoT 1 ^ 2 / (! ~ r)- Sl r- S2 dr - a^ sl a7 2 S2 B(l - si, 1 - s 2 ) . 
Jo 

Since i?(l — Si, 1 — s 2 ) = T(l — Si)T(l — s 2 )/T(2 — Si — s 2 ) we see that the function 
$ is no longer entire. This is to be expected however, because the new polynomial 
f{z\, z%) = (1 + aoiZi)(l + 012.22) has a different Newton polygon, and the new <& 
should contribute to the change of T-factors in the Mellin transform. In fact, when 
«02 = o-ii = we have the formula 

M 1/f ( s ) = $( s )r( S i)r( S2 )r(2 - Sl - s 2 ) = a^a^r^r^ra - Sl )r(i - s 2 ) . 

It is not always the case that all the polar hyperplanes of the gamma functions 
in the representation ([9j) are actual singularities for the Mellin transform M\jf. It 
may happen that the entire function $ has zeros that cancel out some of the poles. 
A very simple example of this phenomenon is provided by the function f(z) = l+z m 
with m > 2. In this case the substitution z m — w leads to the formula 

M 1/f (s) = - / — = — r {s m r 1 - sm) , 

m In 1 + w w m 
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so the polar locus is just mZ. In fact, the entire function <E> from ^ is given by 

, . 1 Tts/rii) r(l — s/m) 

*(*) = rfr^^T ~T-> 

m 1 (s) 1 (to — sj 

and it has plenty of integer zeros. A slight generalization of this example is provided 
by the following result. 

Proposition 2. Let f(z) = 1 + z ai + . . . + z a ", for some linearly independent 
vectors eti, ■ ■ . , ct n € Z™, and denote by S the non-zero determinant det(ajfc). The 
Mellin transform of the rational function 1 // is then given by 

M 1/f (s) = i r«ft,a)) ■ ■ • T((/3 n , b))T(1 - {ft,*)-... - (j3 n , s}) , 

where the tik denote the column vectors of the inverse matrix (a,-fe) . 

PROOF. We make the monomial change of variables z" jl ■ ■ ■ Zn 3n = Wj, so that 
zj = u>f jl • • • Wn 3n and dz/ z — S~~ 1 dw/w. The Mellin transform can then be written 



M 1/f (s) 



dz 1 f W { f US) ---W { n 0n ' S) dw 



n 1 + Z" 1 + . . . + Z Q " Z 8 7 R ™ 1 + Wi + . . . + W n W 



and the latter integral is of a similar form as the one in Example 1. □ 



We point out that the Newton polytope Af of the polynomial in Proposition's 
a simplex with one vertex at the origin, and that its normal vectors fXi, . . . , /i n +i are 
integer multiples of the rational vectors . . . , j3 n and — + . . . + f3 n )- Moreover, 
one has v\ = . . . = v n = and v n +\ = 1. In this case the entire function <f> 
occurring in ^ is therefore of the form 

1 r((/?i,s)) T(({3 n ,s))r(l-(f3 u s)-...-(p n ,s)) 



<5r ((^l) s )) r ((M«, s >) r(l + (Hn+1,8)) 



5. Mellin transforms and coamoebas 

Let us return for a moment to the one-variable Mellin transform 
„ , , f°° z s dz [°° e sx , 

where we assume, as before, that the polynomial / does not vanish on the positive 
real axis and that the real part of s lies in the interior of the Newton interval A f . 
Our first claim is now that the value of the above integral remains unchanged if the 
set of integration is rotated slightly. In other words, for \9\ small enough one has 
the identity 

r°° z s dz _ r z s dz _ r°° e < x+ie ^> 

J J{z)~~ Ja^-Hb) W) ~ ~ J-oo f(e x+id ) 

To verify this, we perform an integration along a closed path starting at the origin, 
then running along the positive real axis to the point R, continuing along the circle 
\z\ = R to the point Re 10 , and then going straight back to the origin, see Figure 4 
below. Since 9 is close to zero, the denominator / has no zeros in the closed sector 
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with arguments between and 9. By the residue theorem the integral over the 
closed contour is therefore equal to zero, and since the integrand decreases fast 
when \z\ — > oo, the integral over the circular arc Cr can be made arbitrarily small 
by choosing R large enough. The integrals along the two infinite rays are thus equal 
as claimed. 




Figure 4. The contour of integration in the residue computation. 



From the above argument we see that the directional Mellin transform coincides 
with the standard one as long as the two directions 9 and belong to the same 
connected component of the coamoeba complement M.\A' f . Furthermore, it is clear 
that the Mellin integral over Arg -1 (#) converges for every choice of 9 outside the 
coamoeba A!j. A similar residue computation as above then again shows that the 
directional Mellin transform only depends on which connected component of K\«4^ 
it is that contains 9. 

Turning to the general case n > 1, there are two important differences to be 
observed. On the one hand we recall the condition in Theorem[T]that the polynomial 
/ should be completely non- vanishing on R™ = Arg _1 (0) in order for the integral 
to converge, and on the other hand we note that the coamoeba A'f is in general 
not a closed set. The following result connects these two facts, and it allows us to 
define the directional Mellin transform 

z s dz _ f e ( a > x+ie '> 

for each argument 9 that does not belong to the closure A'f. 

Theorem 3. For any 9 G K." \ A'j the polynomial f is completely non-vanishing 
on the set Arg _1 (0). 

PROOF. For any given argument vector 9 we can consider the new polynomial 
f e (z) = f( e l9l zi, . . . , e i9n z n ). Observe that A' fB + 9 = A' f , so that G A' fB if and 

only if 9 € A'f, and also that f e is completely non- vanishing on Arg _1 (0) if and 
only if / is completely non- vanishing on Arg _1 (6'). This means that it actually 
suffices to prove the theorem for the special case 9 = 0. 
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Assume then that / is not completely non- vanishing on the set Arg 1 (0), so 
that for some face T one has G A' r . In other words, there is an Xo G M" such 
that /r(e x °) = 0. We must show that belongs to the closure A'j. This is obvious 
if r = A/, so we can assume that dimT < n — 1. Choose a vector fi £ Z n and an 
integer v g Z such that (/x, a) = v for a e T and (//, a) < i> for a € A/ \ T. Writing 
9r = f ~ fr we then have 

aer aer 

and 

aEAj\r aeAy.\r 

where & Q = a Q e^ 0:Q ^ and c a = (fx, a) — v > 0. Now let e > be given. Choose a disk 
D E of radius e centered at so and contained in a complex line on which the function 
w i y fr(e w ) does not vanish identically. Then translate this disk along the real 
space, so that D £ — tfj, is a disk centered at the point Xg — i/i for some large positive 
number t. Since fr(e w ) is non-zero on the boundary of D e we have l/rCe 1 ")! > 6 > 
for w € 9£> £ . This means that l/rCe™)! > 8 e~ tv on the translated circle dD e — i/i. 
Taking t large enough, we also have \gr{e w )/ fr(e w )\ < 1 on dD £ — tfi, that is, 
\gr{e w )\ < \fr{e w )\. Rouche's theorem then tells us that f(e w ) = f T (e w ) + g r {e w ) 
has a zero w e in the disk D E — t/j,. So z E — e We belongs to the hypersurface f(z) = 0. 
But we also know that |Arg(z £ )| = |Imui e | < e, and since e was chosen arbitrary 
we conclude that € A'*. □ 

Remark. In the above proof we showed that all the facial coamoebas A' T are 
contained in the closure A'j of the main coamoeba. It is a fact, proved by Johansson 
[12j and independently by Nisse and Sottile |14j , that one actually has an equality 



u ^ 



A' 



rcA, 



Using Theorems [T] and |3] we can now define a directional Mellin transform (18) 
for any 8 in the complement K™ \ A'^. Just as in the one- variable case discussed 
earlier in the section, the various Mellin transforms will in fact be equal for all 
9 that belong to the same connected component of W l \ A'j. This can be seen 
by connecting two different values of 9 through a polygonal path such that along 
each edge of the path only one component 9^ is being changed. The invariance of 
the Mellin transform under such a move is then a consequence of the one-variable 
argument. 

In order to put our next theorem in a proper perspective, it seems appropriate 
at this juncture to recall some known facts about amoebas and Laments series 
of rational functions. A reference for these results is [6]. Associated with each 
connected component E of the amoeba complement K ra \ A/ is a Laurent series 
representation 

of the rational function 1//. The coefficients of the series are given by the integrals 
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7 a dr f Ja,x+i6) 



(2™)™ Aog-Hx) /(«) 2 ./[-*,*]» /(e x+ie , 

where x is any point in the connected component E. Each such Laurent series will 
converge in the corresponding Reinhardt domain Log^ 1 (i?). We stress the fact that 
the amoeba Aj is always a closed set, so in contrast to the case of coamoebas, there 
is no need to take the closure of the amoeba. 

The following result about coamoebas and Mellin transforms provides a prac- 
tically perfect analogy to the above picture for amoebas and Laurent coefficients. 

Theorem 4. For any connected component E of the coamoeba complement M n \A'f 
there is an integral representation 

i r „ B 



(19) I _M? /f (s)z- s ds, 



which converges for all z in the domain Arg (E) . Here a is an arbitrary point in 
int A f and 

If z s dz 1 f e ( s > x + ie ) 

(20) M 1/f (s) = —J^_^ —- = j^J^ JT-^r^i 

with 9 being an arbitrary point in the component E. 

PROOF. From Theorem[3]and (an obvious generalization of) Theorem[l]we see 
that the integral (20) converges, and from the discussion preceding Theorem [3] we 



also know that the value of (|20|) is independent of the particular choice of point 
G E. 



In order to prove the identity ( 19 1 it suffices to verify that, for all s = a + it 
such that a G int(Ay), the function x — > e ( s ' x+ie ^ / f(e x+l6 ) is in the Schwartz space 
6>(IR n ) of rapidly decreasing functions. Then the result follows from well known 
facts about inversion of Fourier transforms, see Thm 7.1.5 in |10j . 

For simplicity, and without loss of generality, we assume that 8 = 0. We have 
\e^ s ' x "> I f(e x )\ = e^''^ /\f(e x )\, and from the inequality which we established 
in the proof of Theoremlll we see that e^'^ / f(e x ) is an exponentially dercreasing 
function. It remains to verify that all its partial derivatives have the same property. 
Computing a typical derivative, we get 

d (e^^\ a k e^' x ^ e< ff+efe > x > f' k (e x ) 



(21) dxXf(e x )) 



Ck^f(e x )J f(e x ) f{e x f 

where f' k denotes the derivative of the polynomial / with respect to z%. Here the 
first term one the right hand side is just a constant times the original function, and 
the second term is of the form 

a k a a e ^ +a ' x ) 



E 



f{e x ) 2 

The Newton polytope of the denominator is A^ 2 = 2 A f, so a+a € int A p for every 
a E A, and hence each term in the sum satisfies the conditions of Theorem [T] This 
means that the derivative ( |21[ ) is a finite sum of functions to which we can apply 
Theorem [T] and the inequality (|6| . By induction this implies that all derivatives of 
e ( s , x ) j f{e x ) decrease exponentially. □ 
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Remark. It is clear that if E' is a connected component oiW n \A!f that is obtained 
by just translating another component E by 27refc, then the corresponding two 
Mellin transforms are related by the simple formula 

Mf/ / ( S ) = e 2 ^Mf // ( S ). 

In general however, the relations between the various Mellin transforms associated 
with different connected components are rather complicated. Furthermore, it is 
worth mentioning that all the connected components of W 1 \ A'j are convex sets. 
This fact follows for instance from the Bochner tube theorem, see [4]. Finally, we 
point out that Theorem [4] can also be proved by using results from Antipova [2] . 

6. Hypergeometry 

In this final section we shall consider the dependence of the Mellin transform, 
and in particular of the entire function $, on the coefficients a = {a a } of the 
polynomial /. In order to emphasize this dependence we are here going to write 
$(a, s) rather than just &(s). The crucial observation will be that, with respect 
to the variables a, the function $ is an A-hypergeometric function in the sense 
of Gelfand, Kapranov and Zelevinsky. More precisely, a i-> $(a, s) satisfies the A- 
hypergeometric system of partial differential equations with homogeneity parameter 
f3 = (-1, -s\, -s 2 , . . . , -s n ). 

Let us recall the structure of the ^4-hypergeometric system. Our starting point 
is the subset A C Z" of exponent vectors occurring in the expression ^ for the 
polynomial /. We introduce a numbering an, ... , a at of the elements of A, with 
each ak = (ctifc, • ■ ■ , a n k) € Z™. Abusing the notation slightly, we write A also for 
the (1 + n) x TV-matrix whose column vectors are (l,a}~). For any vector v G Z n 
we denote by v + and v~ the vectors obtained from v by replacing each component 
Vk by max(v/ c ,0) and max(— Ufc,0) respectively, so that v = v + — v~ . 

Definition 2. Let A denote a subset {a±, . . . ,q:n} C Z n and the associated (I + 
n) x N -matrix as above. The A-hypergeometric system of differential equations 
with homogeneity parameter f3 £ C" is then given by 

□ b F(a)=0, be1 N ,Ab = 0, and E^F(a)^0, j = 0,l,...,n, 
where the differential operators Ob and Ej are given by 

( d \b + ( d \b~ , g A d 

fc=i 

An analytic function F that solves the system is called A-hypergeometric with ho- 
mogeneity parameter f3. 

Remark. We are assuming N > 1 + n, and as soon as this inequality is strict there 
are of course infinitely many vectors b satisfying Ab = 0, but it is a known fact, see 
[15] . that the system is in fact determined by a finite number of operators 

Let us now, for a given choice of coefficients a, consider an entire function 
s i — y $(a, s) as described in Theorems [2]and |4j We want to study what happens 
when we start varying a. Recall from [J] and [5] the notion of the principal A- 
dctcrminant Ea, also known as the full A-discriminant. It is a polynomial in the 
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variables a, with the property that its zero set C C N contains the singular locus 
of all A-hypergeometric functions. 

Theorem 5. Take a 6 C N \ J^a and let E be a connected component o/R n \ A.'*, 
with f being the polynomial f(z) = a\z ai + . . . + a^z aN . Also take s G C ra with 
Re s £ int A t . Then the analytic germ 

$(o,a) = r r/ ^ — — , 9eE, 

/ias a (multivalued) analytic continuation to (C N \ Ha) x C n which is everywhere 
A-hypergeometric in a with varying homogeneity parameter (—1, — s±, . . . , — s n ). 

PROOF. First of all it is clear that 8 will be disjoint from A!f also for poly- 
nomials / with coefficients a& near the original ones, say in a small ball B(a), so that 
the integral does indeed define an analytic germ $(a, s). From (a straightforward 
generalization of) our Theorem [5] we also know that $ is extendable as an entire 
function with respect to the variables s. In other words, we already have an analytic 
extension of $ to the infinite cylinder B{a) x C n . 

Let us next verify that $ is an A-hypergeometric function with the correct 
homogeneity parameter. When doing this we first fix s at an arbitrary value with 
Re s £ int A f , hence in particular away from the polar hypcrplancs of the gamma 
functions. Then the function in front of the integral is just a non-zero constant and 
we can deal directly with the integral, by differentiation under the integral sign. 

Notice that the condition that Ab = amounts to the two identities = \b~\ 
and (b + ,a) = (b~,a), were we have used the shorthand notation |6 ± | = J2 an< ^ 

(b ± ,a) = Computing iterated derivatives of the integrand 1// in the 

Mellin integral we get 

/ f) \ b ± 1 r ( b± ' a ) 

(L) 1 = (-ni b ± i|6±ii - 

\daJ J2 a kz ak (£akZ a ") 1+ \ b± \ ' 

and since here the right hand side is independent of the choice of sign in b ± , SO IS 
the left hand side. This means that = 0, and hence we also have = 0. 

It is obvious that <& is homogeneous of degree —I with respect to the variables 
ak ■ To check the other homogeneities one can integrate by parts in the integral. As 
in our proof of Theorem [2] this can be efficiently done by dilating the variables by 
means of a parameter A. For example, making the dilation Zj > Xzj we get 



/ 

JAi 



z s dz , „ f z s dz 

= A 



I Arg- 1 (6) J2k a kZ ak Z Arg-He) Efc X* ik a k Z a * Z 

Differentiating both sides of this identity with respect to A and then putting A = 1, 
we find that 

N q 
k— 1 

and hence Ej$(a, s) = 0, with f3j — —sj as claimed. 

We have thus established that $ is an A-hypergeometric analytic function in the 
product domain B(a) x (int Af + i R n ), and by uniqueness of analytic continuation 
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its extension to the cylinder B(a) x C" will remain A-hypergeometric. Next, by the 
general theory of A-hypergeometric functions one has, for each fixed s, a (typically 
multivalued) analytic continuation of o 4 $(a, s) from B(a) to all of \ XU- 
Well known results on analytic functions of several variables then tell us that these 
continuations will still depend analytically on s, so we have achieved the desired 
analytic continuation to the full product domain (C N \ E^) x C™. The uniqueness 
of analytic continuation again guarantees that $ will everywhere satisfy the A- 
hypergeometric system with the homogeneity parameter (— 1, — s\. . . . , — s n ). □ 

Related integral representations of A-hypergeometric functions have been con- 
sidered by several authors, see for instance [9] and [3]. It is probably instructive to 
examine a concrete special instance of the above theorem, and we choose to present 
the case of the classical Gauss hypergeometric function. 

Example 3. Take A — {(0, 0), (1, 0), (0, 1), (1, 1)} to consist of the four corners 
of the unit square in the first quadrant. It is easy to check that in this case £a 
is given by the equation Ea(cl) = aia 2 a^ai{a\ai — 0203) = 0. Then consider the 
polynomial f(z) = a\ + 02^1 + 02^2 + a4.z\z 2 together with its associated Mellin 
transform 

z-i 1 Zn 2 dz\dz 2 



(22) M 1/f (a,s) , , 

Jo Jo ai + a 2 zx + a 3 z 2 + a4Zxz 2 z x z 2 

Let us compute this transform for simplicity first in the case a\ = a 2 — = 1. 
Writing / as (1 + z 2 ) + (1 + a^z^zi, we can use formula (16) and first perform the 
integration with respect to z\ . This yields the expression 

/•oo 

r(ai)r(l-«i) / (l + z 2 ) Sl - 1 (l + a 4 z 2 )- Sl z s 2 2 ~ 1 dz 2 
Jo 



for the Mellin transform ( 22 ) . Re- writing the integrand and expanding in a power 
rat the above integral equals 

1)^ -1 z s r 1 dz 2 = ^ ^r(l-sQ (fl4 _ 1)fc [°° z? +k ~\ 
so using the formula 



series we find that the above integral equals 

, («4 ' ■ - \ " - - _ \ • .1. ' ( 

1 + ^2 ) i + z 2 ^ o r(i-, Sl -fc)fc! 1 ^' Jo (i + ^) 1+fe: 



w* dw 1 , , 

n , u+fc — = t, r t r i + k - 1) 

(1 + w) L+h w fc! 



we find that (22) can be expressed as r(si)T(l — Si)r(s 2 )r(l — s 2 ) times 



\- T(l- Sl )T(s 2 + k) r( Sl + k)T(s 2 + k) k 

^ o r(i- Sl -k)T( S2 )(k\y (a4 L) ^ r( Sl )r( S2 )(fc!) 2 [L ° 4j ' 

In other words, we have shown that $(1, 1, 1, 04, s±, s 2 ) = 2 Fi(si, s 2 ; 1; 1 — a 4 ). 
Using the homogeneities if $, corresponding to the row vectors of the matrix 



A : 
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and the homogeneity parameter (— 1, — s\, — s 2 ), we then easily recover the more 
general formula 

(23) 1>(a, a) = a^+^a^a^ 2 F l ( Sl ,s 2] 1; 1 - — ) . 

0203 

Notice that the three singular points 0, 1 and 00 for the Gauss function correspond 
to the factors 0104 — 0.203, 0104 and 0203 of the principal A-determinant Ea- 

The A-hypergeometric system consists in this case of the single binomial equation 

8 2 8 2 \ 

together with the three homogeneity equations 

{(a\d/dai + a 2 8/8a 2 + 038/803 + 048/804 + 1) <&(a, s) = 0, 
(a 2 8/8a 2 + 048/804 + s x ) <f>(a, s) = 0, 
(0,38/80,3 + 0,48/80,4 + s 2 )<5>(a,s) = 0. 

Let us end by considering what happens as one of the variables Ok vanishes. 
From the Gauss hypergeometric theorem one knows that, for ax + a 2 < 1, there is 
an identity 2 Fx(sx, s 2 ; 1; 1) = T(l — Si — s 2 )/(r(l — si)r(l — s 2 )), so setting 04 = 
in the above formula ( p3| we get 

81+^2-1 „-»i „-* 3 T(l - s x - 33) 



<J>(ai, a 2 , a 3 , 0, sx, s 2 ) = a* 1 82 a 2 Sl a 



r(i- Sl )r(i-s 2 ) ' 

This of course fits beautifully with the fact that the Mellin transform M\ 1 f for the 
polynomial f(z\, z 2 ) = ai + a 2 z\ + 0.32:2 is equal to 



S1+S2 — 1 — Sx — S2 



Oi a n a 



- S2 r( Sl )r( S2 )r(i- Sl - S2 ), 



compare with formula ( 16 1 above. 

To treat the case a\ — we can use Euler's hypergeometric tran sfor mation 
2-Fi(a, 6; 1; z) = (1 — z) x ~ a 2 Fx{^ — a, 1 — b; 1; z) and re-write formula (23) as 

<I>(a, s) = a 2 2 a 4 2 2 ^i(l - si, 1 - s 2 ; 1; 1 J. 

02^3 

Again using the hypergeometric theorem of Gauss, now for (1 — <j\) + (1 — a 2 ) < 1 
or equivalently <j\ + a 2 > 1 , we find that 

*2-i„si-i r ( s i + s 2 - 1) 



$(0, a 2 , a 3 , a 4 , Si, s 2 ) = a s 2 2 a^ 1 a 4 



r(*i)r(« a ) ' 

which can be seen to concord with the formula for the Mellin transform M\if with 
f{z\, z 2 ) — a 2 z\ + a3Z 2 + a4Z\z 2 . Finally, the other two Euler transformations 

2 Fx(a,b;l;z) = (l-z)- a 2 Fx(a, 1-6; 1; z/{z-l)) = {l-z)- b 2 Fx(l-a, 6; 1; z/(z-l)) 

similarly lead to the formulas 

$(ai,0,a 3 ,a 4 ,si,s 2 ) = a s t 2 ~ 1 a^' 82 a 

and 

$(ai,a 2 ,0,a4,si,s 2 ) = a\ 



_„-. i^- S2q - s1 r( S2 -si) 
3 fl4 r(i- Sl )r( S2 )' 



2 4 r( Sl )r(i- S2 ) 
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